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ABSTRACT 

"qq". Tides raised on a planet by its host star's gravity can reduce a planet's or- 

O , bital semi-major axis and eccentricity. This effect is only relevant for planets 

Psj . orbiting very close to their host stars. The habitable zones of low-mass stars 

are also close-in and tides can alter the orbits of planets in these locations. We 

calculate the tidal evolution of hypothetical terrestrial planets around low-mass 

stars and show that tides can evolve planets past the inner edge of the habit- 

. able zone, sometimes in less than 1 billion years. This migration requires large 

p ,! eccentricities (> 0.5) and low-mass stars ( <O.35M0). Such migration may have 

Q ! important implications for the evolution of the atmosphere, internal heating and 

^— > ! the Gaia hypothesis. Similarly, a planet detected interior to the habitable zone 

ci ', could have been habitable in the past. We consider the past habitability of the 

recently-discovered, ~ 5 M^ planet, Gliese 581 c. We find that it could have been 

habitable for reasonable choices of orbital and physical properties as recently as 

O ! 2 Gyr ago. However, when we include constraints derived from the additional 

\Q . companions, we see that most parameter choices that predict past habitability 

O I require the two inner planets of the system to have crossed their mutual 3:1 

L^ ! mean motion resonance. As this crossing would likely have resulted in resonance 

OO ! capture, which is not observed, we conclude that Gl 581 c was probably never 

habitable. 
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1. Introduction 

After the discovery of extra-solar planets, the detection of a habitable, terrestrial (rocky) 
planet became a real possibility. Of particular interest are planets in orbit about low-mass 
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stars ( <O.5M0). The lower stellar mass makes detection of lower mass planets easier and 
leads to habitable orbits with smaller semi-major axes (which also increase the likelihood 
of detection by current technologies). The recent detection of putative 5 and 8 M^ planets 
around an M star (stellar mass M^ = O.SIMq) (Udry et al. 2007) near the circumstellar 
"habitable zone" (HZ) (Kasting et al. 1993; Franck et al. 2000; von Bloh et al. 2007; Selsis 
et a/. 2007; hereafter S07) suggests that Doppler surveys {i.e. Butler et a/. 2006; Naef et 
al. 2007) are becoming sensitive enough to detect terrestrial planets in the HZ of low-mass 
stars. Not surprisingly, M stars now appear to be attractive targets for detecting habitable 
planets (Tarter et a/. 2007). 

The habitable zone is defined, following Kasting et al. (1993), as the range of distances 
from a star at which liquid water could be stable at the surface of a planet. The definition 
usually implicitly assumes that the star is the main source of heat. Thus, it excludes locations 
where other sources of heat might be available, such as moons, like Jupiter's Europa, where 
tidal heating is a major factor maintaining near-surface liquid water. In that sense, the 
phrase habitable zone, as usually used in astrophysics, is a misnomer. Nevertheless, here we 
use the conventional definition, while explicitly emphasizing that life could exist outside this 
zone. 

The detection of Gl 581 c and d (Udry et al. 2007) has demonstrated that planets can 
form near the HZ of low-mass stars. For such stars the habitable zone is so close to the star 
that tidal effects may be important (Mardling & Lin 2004). Tides raised between the star 
and planet introduce torques on the planet resulting in orbital evolution, usually decay of 
semi-major axis, a, and eccentricity, e. Here we quantify the effects of tidal forces in such 
zones. We find that Earth-sized planets with large initial eccentricities could migrate out of 
the HZ of low-mass stars in less than 1 billion years. 

The proper definition of an HZ is not obvious. As described in S07, the boundaries of 
the HZ are a complex function of both stellar and planetary properties. For any individual 
planet, the definitions we choose may not be appropriate. Nonetheless, we will show a general 
trend that must be considered when assessing a planet's habitability, both in the past and 
in the future. 

We consider the evolution of terrestrial planets beginning from a location that is almost 
assuredly habitable to a region that is probably not habitable. We call the time to complete 
this traversal the "habitable lifetime", thl- For relatively massive stars ( >0.35M7)), the HZ 
is far enough away from the star that tides cannot cause significant tidal evolutioio. However 
for planets with relatively large eccentricities ( >0.5) orbiting low-mass stars ( <O.35M0) 
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tides can pull a planet out of the HZ in less than 1 Gyr. The large number of low-mass stars 
and significant fraction of known exoplanets with large e suggest that tides may limit the 
total number of habitable planets in the galaxy. 

We assume that potentially habitable planets can form with a range of masses and 
orbits in the HZ of low- mass stars (see Chambers [2004]; Raymond, Quinn & Lunine [2007] or 
Lissauer [2007] for a discussion of the formation of such planets) . Of course, the formation of a 
planet in the HZ does not mean a planet is, will be, or has been habitable. Many other factors 
are required such as favorable composition (Lissauer 2007), a thick atmosphere, substantial 
mass for sustained plate tectonics (m >0.3M®) (Williams, Kasting & Wade 1997; Raymond, 
Scalo & Meadows 2007), which on Earth drives the carbonate-silicate cycle (Walker, Hays & 
Kasting 1981), and substantial water content {e.g. Raymond, Quinn & Lunine 2004), both 
for life and to facilitate plate tectonics (Regenauer-Lieb, Yuen & Barlund 2001). Nonetheless 
we will call planets "habitable" if their orbital properties are such that they permit water to 
be stable on the surface. 

The recently announced Gl 581 system contains a planet, c, that is probably just interior 
to the HZ (S07). In Table 1 we list the properties of the planets in this system; rrip is the 
planet mass, P is the orbital period, a is the semi-major axis, e is the eccentricity, w is the 
longitude of periastron, and Tperi is the time of periastron passage in Julian Days. Tidal 
theory predicts planet c orbited with larger values of semi-major axis and eccentricity in 
the past. Therefore it may have been habitable in the past, but tides subsequently moved 
the planet into an uninhabitable orbit. We find that if planet c was the only planet in the 
system, plausible physical properties predict it was habitable. However, when we consider 
constraints from the mutual interactions of the additional planets, we find that planet c has 
likely never been habitable. 

Tides strong enough to induce significant orbital evolution can also heat the interior 
of a planet through tidal working. The same type of heating causes Jo's intense vulcanism 
(Peale, Cassen & Reynolds 1979). Therefore we might expect an increased heating rate 
on the planets we consider here as well. The orbital evolution we consider could also have 
implications for the evolution of a habitable planet's atmosphere. In some cases, the orbit- 
averaged flux can increase by a factor of 2 as the orbit shrinks (assuming constant stellar 
luminosity) . 

In § 2 we review tidal theory, present our deflnition of the HZ, and describe how we 
determine the habitable lifetime. In § 3 we present our results for hypothetical terrestrial 
planets. In § 4 we describe how our model applies to the Gl 581 c planet. Finally, in § 5 we 
draw our conclusions and discuss directions for future research. 
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2. Methodology 

2.1. Tidal Theory 

We evaluate the effects of tides using the classical, second order equations assembled by 
Goldreich & Soter (1966) (see also Kaula 1964; Greenberg 1977): 

da ^ _ / ^ v^G^^^2 9y^G/m^m^~.^_,,/, 

dt \ rripQp ^ 2 Q', / ^ ' 



dt V 2 mpQp 16 Q'^ 

where a is semi-major axis, G is Newton's gravitational constant, m^ is the mass of the 
planet, k is the planet's Love number, Qp is the planet's tidal dissipation function, Q'^ is 
the star's tidal dissipation function divided by two-thirds its love number, Rp is the planet's 
radius, and /2* is the star's radius. The first terms in these equations are due to the tide 
raised on the planet by the star, and the second terms are due to tides raised on the star 
by the planet. Eqs. ([I]-[2]) are only valid if the planet's orbital period is less than the star's 
rotation period, otherwise the signs of the second two terms will change. For the cases we 
consider, this condition is irrelevant because the second terms of Eqs. ([T]-[2]) are 6-7 orders 
of magnitude smaller than the first. This difference is largely due to the mass difference 
between a terrestrial planet and a star. Note that tidal theory predicts decay for both a and 
e. As tides slowly change a planet's orbit, the planet may move into or out of the HZ. 

The equations above were derived under the implicit assumption that e values are small. 
Here we apply them to cases that include conditions with large e values. Thus higher-order 
corrections may potentially be important. Several models of tidal evolution that are not 
restricted to small values of e have been derived (Hut 1981; Eggleton, Kiseleva & Hut 1998; 
Mardling & Lin 2002). However, such models must include specific assumptions about how 
a body responds to the ever-changing tidal potential. The tidal forcing involves components 
with a various frequencies, but the nature of the response of a planet or star is uncertain. 
It is not clear how the tidal phase lag depends on the frequency of each harmonic. This 
issue becomes increasingly critical when e is large and a wide range of different comparable- 
amplitude frequencies are involved. For a model valid to a high order in e, the validity of 
parameterizing tidal effects through a single constant Q is questionable. Moreover, tidal 
heating could lead to core/mantle melting and a time- dependent Q value. Not enough 
is known about the actual response of real bodies to evaluate these higher-order effects. 
Therefore, to the best of current knowledge, Eqs. (1 - 2) above are reasonable representations 
of tidal evolution. As the theory matures, the effect of higher-order terms could be included 
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and this study should be revisited. Note that the higher order theories that have been derived 
all predict faster evolution than Eqs. (1 - 2) (all the higher order terms have the same sign 
as those in Eqs. (1 - 2), so their predictions for the habitable lifetime, thl^ would be even 
shorter than those we derive below. In that sense, our approach provides a conservative 
estimate of the tidal effects on terrestrial planets. 

Eqs. ([I]-[2]) are only valid when the star's orbit-averaged torque on the planet's tidal 
bulge is zero. When this situation occurs, the planet is said to be "tidally locked" to its 
star. For a planet on an eccentric orbit, the torque is greatest at pericenter at which point 
the star is moving with an angular velocity greater than the mean motion. Therefore the 
rotational frequency of the planet Vt is larger than its orbital frequency n (tidal locking does 
not mean that one side of the planet faces the star for all times, this situation only occurs 
for a circular orbit). To second order in e these two frequencies are related by the following 
equation, 

^ = n{l + -e') (3) 

(Murray & Dermott 1999; Goldreich 1966). Note that this equation neglects effects of 
a permanent quadrupole moment (we ignore this effect as their is currently no hope of 
measuring it in an exoplanet). Even for modest eccentricities, the spin and orbital frequencies 
can be significantly different due to the coefficient in front of the e^ term. A planet with an 
eccentricity of 0.32 has a spin period half that of its orbital period. 



The time for a planet to become tidally locked is given by 

AbGM^kr. 






p 



(Goldreich & Soter 1966; Rasio et al. 1996), where Af2 is the difference between the initial 
spin, and the spin given by Eq. ^. For an Earth-mass, terrestrial planet {k = 0.3, Qp = 100) 
planet orbiting a 0.3 Mq star with e = a = 0.1 AU, this timescale is about 10^ years. We 
consider planets initially with much larger eccentricities, but eccentricity probably does not 
affect this timescale significantly (Peale 1977). This time is very short compared to other 
timescales associated with planet formation, and we may assume the planets are tidally 
locked to the star as soon as they are formed. 

In order to determine thl for possible planets, we consider a range of values of M^, nip, 
a and e. We will set Q'^ = 10^'^, Qp = 21.5, and k = 0.3, reasonable choices based on limited 
observational constraints for stars (Ogilvie & Lin 2007; Jackson, Greenberg & Barnes 2007) 
and the Earth (Dickey et al. 1994; Mardling & Lin 2004). However, the time-averaged value 
of Q of the Earth is probably closer to 100 (Lambeck 1977), which is similar to the values of 
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other terrestrial planets in the Solar System (Yoder 1995). Nonetheless we choose Qp = 21.5 
to be consistent with past work (Mardling & Lin 2004). 

From our choices of M^ and nip we must determine radii. For the planet, we will use 

the scaling relationship 

7? / m \ 0.27 

(Valencia, O'Connell & Sasselov 2006). For the stellar radius, we use an empirical relation 
derived from observations of eclipsing binaries (Gorda & Svechnikov 1999): 

logio^ = 1.03 logiov7^ + 0.1, (6) 

where Rq and Mq are the solar radius and mass, respectively. Note that Eq. (Q is only 
valid when M* <Mq. 



2.2. The Eccentric Habitable Zone 

For circular orbits, the HZ is defined as the range of orbits in which a terrestrial planet, 
with an atmosphere favorable for habitability, can support liquid water on the surface (Kast- 
ing et al. 1993; S07). However, extra-solar planet orbits are often highly eccentric. For these 
orbits, Williams & Pollard (2002) showed that, despite large variations in surface tempera- 
ture, the key for long-term chmate stabihty is the orbit-averaged flux, 

(Williams & Pollard 2002; Adams & Laughlin 2006), where F is the orbit-averaged flux and L 
is the stellar luminosity. Note, however, that potential water loss at periastron and potential 
cloud freeze-out at apoastron have not been accounted for, despite their implications for 
habitability. We define the eccentric habitable zone (EHZ) as the range of orbits for which a 
planet receives as much flux over 1 orbit as a planet on a circular orbit in the "classical" HZ 
(Kasting et al. 1993; S07). We must also scale the location of the EHZ to the luminosity of 
the stars we are considering. Lower mass stars have lower luminosities, and hence the EHZ's 
are closer to the star. We therefore define the inner and outer edges of the EHZ, /j„ and lout 
(in AU), respectively, as 

T 1 /9 

kn = {kne - a^nT, ' hrJ^) (—) (1 - e^)-^'\ (8) 

and 



lout = {loutG - cioutT* - boutT^) il—) (1 - ^) ^^^' (9) 



-7- 



where Im and lout are the inner and outer edges of the HZ, respectively, in AU, ImQ and 
loute ^^6 the inner and outer edges of the HZ in the solar system, respectively, in AU, 
ttin = 2.7619 X 10-5 AU/K, k^ = 3.8095 x lO'^ AU/K^, a^ut = 1-3786 x 10"^ AU/K, and 
bout = 1.4286 X 10^^ AU/K^ are empirically determined constants, T^ = T^ff — 5700 K, 
and Lq is the solar luminosity (S07). Terms involving T* relate stellar intensity at different 
wavelengths and atmospheric windows of predicted habitable planets. Eq. ([8] - [9]) are only 
valid for T^jf < 3700 K. T^fj is the effective temperature of the star, given by 

where a is the Stefan-Boltzman const ant. Therefore a planet with a given e is on the inner 
edge of the EHZ if a = /,„. Note that as e changes so do the locations of /„ and lout- 

The values of /^q and loute ^^^ therefore the key parameters in the identification of the 
edges of the EHZ. We consider three criteria identified in S07: 1) 0% cloud cover, 2) 50% 
cloud cover, and 3) 100% cloud covecl. S07 give values of I^q of ~ 0.89, ~ 0.72, and ~ 0.49 
AU for the three possibilities, respectively. For the outer edge, S07 give loute = 1-67, 1.95, 
and 2.4 AU, for the three cloud cover models, respectively. We arbitrarily choose the 50% 
cloud cover case to be the limits of the EHZ. This choice is the middle of the possibilities, 
and roughly corresponds to the cloud cover on the Earth. 

The stellar mass is related to the luminosity in the following way, 

A = 4.101/i^ + 8.162/^2 + 7.108/i + 0.065, (11) 

where A = logio(-^/-^0) and /i = logio(M*/M0) (Scalo et al. 2007). With these definitions, 
we plot the shape of the EHZ as a function of stellar mass and initial eccentricity, Cq, in Fig. 
[U As both Co and M* increase, the location of the EHZ moves out and grows wider. If, for 
a certain star, Tg// < 3700 K, we set it to 3700 K in Eqs. ([8] -[9]) as done in S07. 



2.3. Determination of the Habitable Lifetime 

The habitable lifetime, thl, depends on the initial orbit of a planet. Rather than explore 
all possible orbits, we will consider one type of starting location: the inner 0% cloud cover 
boundary {i.e. ImQ = 0.89 AU). Although this choice of initial condition is arbitrary, it 
demonstrates how tides will change the orbits, and, hence, habitability. As stated in the 



^Note that S07 also give values of linQ and loutQ based on assumptions regarding the past and present 
habitability of Venus and Mars coupled with stellar evolution, but we will ignore those models here. 



previous section, we consider a planet uninhabitable when it crosses the 50% cloud cover 
boundary {linQ = 0.72 AU). We vary the planet's initial eccentricity Cq and M* in increments 
of 0.01 and 0.01 Mq, respectively, and consider four different planet masses: 5, 1, 0.5, and 

0.1 Me. 

In order to calculate thl we numerically integrate Eqs. (P - E]) and ([8]) with a timestep 
of 10,000 years. Convergence tests showed that this timestep is two orders of magnitude 
smaller than necessary to produce reliable results. We integrate the tidal equation until 
either a < /j„(50% cloud cover) or thl > 10 Gyr. 



3. Habitable Lifetimes 

In this section we present results for thl as defined in § 2.3. Fig. [2] shows contour lines 
of Thl for planets with masses of 5, 1, 0.5 and 0.1 M^. There are two prominent features 
in each panel. First is the abrupt change in thl from 10 Gyr to less than 1 Gyr between 
Co = 0.53 and 0.54. This feature occurs because the time to circularize the orbit is less than 
10 Gyr. When e reaches zero, the evolution of a effectively stops, and the planet is stuck 
on one side of the boundary or the other. Orbits with cq <0.53 become circular with a just 
larger than /j„. However, when cq >0.54, the orbits circularize with a slightly less than /«„. 

These different types of evolution are shown in Fig. [31 Note that the evolution of both 
cases is very similar, but the final relationship between a and /,„ is such that one planet 
remains in the EHZ, but the other does not. Note that the EHZ changes because e is 
changing. The difference is therefore not due to some fundamental feature of tidal theory 
(other than changes in semi-major axis become negligible for small eccentricities), rather our 
definition of the EHZ leads to the different values of thl- 

The apparently sudden transition from long to short habitable lifetimes can be shown 
by considering how a and e are related: 

a \ 2 2 



In - =e^-e, (12) 

Vao/ 

where ao and cq are the initial values of a and e, if the effects of the tide raised on the star 
by the planet is negligible (which is the case here) (Jackson, Greenberg & Barnes 2007). If 
we assume e = {i.e. the tidal evolution has effectively ended, and t = oo), a = Zj„(50% 
cloud cover) and ao = /m(0% cloud cover), then the value of Cq that solves Eq. flT2l) should 
correspond to the critical value (cq ~ 0.535) seen in Fig. [2] (note that ao is also a function 
of Co and Eq. [12] must be solved numerically). In the cases we consider, this equation is 
independent of mass. The values of a and ao have the same dependence on luminosity (and 



-9- 



hence mass, c.f. Eq. [H]), and, for the stars we consider, we have set Tg// = 3700 in order 
to follow the example of S07. Therefore we expect the critical value of Cq to be independent 
of both stellar and planetary mass. Solving Eq. (IT^ for the values of /j„ we have chosen, we 
find the critical value of Cq is 0.53490, in agreement with Fig. [2J 

The second feature of the contours is the downturn that occurs for eo >0.65, in other 
words, thl gradually increases with cq. This dependence follows from the change in the EHZ 
at larger values of Cq, as shown in Fig. [TJ At these values not only is the EHZ wider, but Oq 
is also larger. Therefore planets' orbits will evolve more slowly (note the a dependence in 
Eqs. P-[2]), and have further to go before they reach /j„(50% cloud cover). 

We also see in Fig. [2] that lower mass planets tend to have longer habitable lifetimes. 
This relationship follows from the dependence of da/dt on mp in Eqs. ([I]-[2]). Low- mass stars, 
like the ones considered here, are likely to host low-mass planets (Laughlin, Bodenheimer & 
Adams 2004; Ida & Lin 2005; Raymond, Scalo & Meadows 2007), and therefore disks around 
M stars may preferentially form planets with relatively long habitable lifetimes. However, 
they may not contain enough volatiles for life because the impacts that deliver them may 
have large enough energies that the volatiles will be lost to the planet (Lissauer 2007). 

Models of the evolution of habitable atmospheres have focused on planets with static 
orbits {e.g. Kasting & Catling 2003). Figure [2] shows that such a model may be inappropriate 
for determining habitability and biosignatures of planets around M stars (Segura et al. 2005; 
S07). In Fig. m we plot the orbit-averaged stellar flux for a planet with t^l = 4.5 Gyr 
orbiting a 0.2 Mq star (assuming the stellar luminosity is constant). This figure also shows 
how the spin frequency of the planet Q evolves relative to its mean motion n, see Eq. (|3]). 
The tidally evolving planet does not orbit such that only one side faces the star until the 
orbit has become circular at t ~ 6 Gyr. 

Note as well that there exists a family of orbits with thl ~ 4.5 Gyr, about the age of 
the Earth. Therefore these calculations suggest that if a planet with the correct ingredients 
for life formed with that orbit, and life subsequently developed in a manner similar to the 
Earth, that life would be eliminated by the ultimate global warming event: The passage of 
the planet through the inner edge of the HZ due to tidal evolution. 



4. Application to the Gl 581 System 

In this section we apply the ideas of § 2 to the recently announced Gl 581 planetary 
system (Udry et al. 2007). This system contains three planets with orbits and masses listed 
in Table 1. Planet c is probably interior to the 50% cloud cover HZ, and planet d is probably 
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exterior (S07). The age of the star is uncertain. Its kinematics and metalhcity suggest 
that it is at least a few Gyr old (Bonfils et a/. 2005; Delfosse et al. 1998), but its very low 
X-ray luminosity (Delfosse et al. 1998) indicates it may be as old as 8 - 10 Gyr old (Ribas, 
private communication; Selsis, private communication). We therefore tidally evolve planet 
c backward for 10 Gyr by changing the signs of Eqs. ([I]-[2]) (Jackson, Greenberg & Barnes 
2007). First we consider how Gl 581 c would evolve if it were an isolated planet, then we 
consider how the additional companions constrain its evolution and the planets' physical 
properties. 



4.1. Past Evolution of Gl 581 c 

If Gl 581 c were the only planetary companion to its host star, then its past tidal 
evolution is adequately modeled by Eqs. ([I]-[2]). For this particular system, we have specific 
constraints on physical and orbital properties. The planet's radius Re lies in the range 
1.6 <Rc <2-Re, assuming the observed minimum mass is the actual mass (S07; Valencia, 
O'Connell & Sasselov 2006; Fortney, Marley & Barnes 2007; Valencia, Sasselov & O'Connell 
2007a, 2007b; Sotin, Grasset & Moquet 2007), and its best-fit eccentricity Cc is 0.16, although 
any value between and 0.3 is about equally likely (Udry et al. 2007). For tidal parameters, 
we appeal to values in our Solar System. The current values of k and Q for the Earth are 
about 0.3 and 21.5, respectively, from Lunar Laser Ranging (Dickey et al. 1994; Mardling & 
Lin 2004), but over the lifetime of the Earth-moon system, the Q value is probably close to 
100 (Lambeck 1977), similar to the value of Mars (Yoder 1995). The star's mass, radius, 
luminosity and effective temperature are O.3IM0, O.38i?0, 0.0135Lq and 3200 K, respectively 
(Bonfils et a/. 2005). Planet c is currently inside the 100% cloud cover HZ, but such an 
atmospheric condition is unlikely (S07). We therefore seek to identify physical and orbital 
parameters that allow Gl 581 c to have been habitable in the past, i.e. inside the 50% cloud 
cover HZ. 

First we consider the effects of varying eccentricity. In Fig. O we show how long ago Gl 
581 c would have been inside the 50% cloud cover HZ, as a function of its current eccentricity, 
for 2 different, plausible values of Qc, assuming k = 0.3. At low Cc values, the time is longer 
because tidal evolution is slower with smaller e, while at larger values the time is longer due 
to the strong coupling between a and e, and because the edge of the HZ is farther out, c.f. 
Eqs. ([I]-[2]) and ([H]). Therefore there exists a critical value of Cc that minimizes the time 
since the planet was habitable. In this case, that value is 0.32 in which case the planet would 
have been inside the 50% cloud cover HZ as recently as 3 Gyr ago. 

Next we consider how the planet's radius constrains past evolution. For most cases we 



-11- 



considered, the Gl 581 c planet could have been inside the 50% cloud cover HZ for many 
Gyr, assuming a system age of 10 Gyr. In Fig. [6] we plot the evolution of the EHZ and a^ 
for three different values of Rp, and assuming its current eccentricity is 0.16. All the cases 
plotted assume Qc = 21.5, k = 0.3. The shading represents the 100%, 50%, and 0% cloud 
cover models, and white is outside the HZ, as in Fig. [H Even for the extreme case Vc = 2R(^, 
the planet could never have been in the 0% cloud cover HZ. In this case, it is hard to imagine 
how Qc could be so low, since the planet is most likely a "water world" (Raymond, Quinn 
& Lunine 2004; Leger et al. 2004) and for such a planet dissipation of tidal energy would be 
less efficient {Q would be larger). If the system is 10 Gyr old, this case predicts planet c was 
habitable for 8 Gyr and tides would have sterilized planet c 2 Gyr ago. 



4.2. Possible Interactions Between Planets b and c 

So far we have ignored interactions between planets b and c, but as we run backwards 
in time, we need to remember that their mutual gravitational interactions will cause the 
eccentricities to oscillate on a timescale of ~ 10^ years, and we must consider the possibility 
that as their periods change, they could hit a mean motion resonance. The first significant 
resonance would be the 3:1 ratio of orbital periods. If they ever crossed this resonance, the 
planets would have almost assuredly been captured into the resonance (Lee & Peale 2002). 
Such a scenario would invalidate the conclusions of the previous section because the system 
is not observed to be in resonance today. 

Planet b is also subject to tidal forces so we must account for its tidal evolution. Cur- 
rently its outer 3:1 resonance lies at ~ 0.085 AU, and therefore if it did not experience 
significant evolution, then planet c has orbited inside that distance since the system formed. 
If, however, planet b experienced significant tidal evolution, then both planets could have 
migrated such that the 3:1 was always instantaneously beyond the orbit of planet c. The 
value of Q'j, is unknown for the Neptune-mass planet. If the planet is gaseous [k = 3/2), 
then Qfy ~ 10^ (Jackson, Greenberg & Barnes 2007), but if it is terrestrial, then the value is 
probably closer to 100. In Fig. [7] we plot the evolution of at and e^ for values of Qb between 
lO'^ and 10^. We assume the planet has the same bulk density as Neptune, and therefore 
the radius of 3.74i?^. This choice is somewhat arbitrary, but is consistent with the radius 
of the Neptune-mass planet GJ 436 (Gillon et al. 2007; Deming et al. 2007). The evolution 
depends strongly on Qh, and in order for planet c to have been habitable and avoided the 
3:1 resonance, then Qf, <5 x 10"^. 

In fact, Qb must be less than 4 x 10"^ in order to avoid the 3:1 resonance, as shown 
in Fig. O This figure plots the evolutions of b and c (without mutual interactions, and 
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Qc = 21.5, Tc = 1.8i?©), and the instantaneous 3:1 resonance of planet b. We see that in this 
case planet c just misses the 3:1 resonance. 

To this point, we have not considered the secular interactions between the planets. 
Secular interactions cause the eccentricities and orientations of the orbits to oscillate with 
time (see e.g. Barnes & Greenberg [2006] for a review of secular theory). In Fig.[8]we show the 
secular interactions for the currently observed system, excluding planet d. This simulation 
was performed with HNBOD"Y|j, and includes general relativistic effects. Our integration 
reveals that the current values of the eccentricities are near their extrema. We presume that 
using the observed values provides a reasonable estimate of the tidal evolution, however we 
may be overestimating the tidal effects on planet c, and underestimating the tidal effects 
on planet b. A more sophisticated treatment that combines the secular and tidal evolutions 
was beyond the scope of this investigation (see § 5). 

Tidal theory predicts the eccentricities would be larger in the past. The amplitude of 
the eccentricity oscillations (due to gravitational interactions between planets) scales with 
e (among other quantities), and eccentricity oscillations would therefore have been larger in 
the past. Larger eccentricity values can often lead to dynamical instability {e.g. Barnes & 
Quinn 2004; Barnes & Greenberg 2007). Therefore we also require that the mutual, tidal 
evolutions of planets b and c predict dynamically stable orbits. Such a consideration shows 
that if Qc is small enough to exclude a 3:1 resonance capture {i.e. < 4 x 10^), then the orbits 
of planets b and c would have been unstable 10 Gyr ago. This instability arises from the 
large e values the planets had then (e?, = 0.74, Cc = 0.47). 

Even if we begin the backward integration with the two planet's average eccentricity 
(over the secular cycle), planet c was probably not habitable. Planet c's average eccentricity 
is lower than its observed value, and it has probably experienced less tidal evolution than 
predicted above. Therefore it is unlikely it could have ever been inside the 50% cloud cover 
EHZ. On the other hand, b's average eccentricity is larger than the current value, and its 
tidal evolution would have probably been more significant. In that case, tidal evolution is 
more likely to predict dynamical instability. We conclude that the requirements that planets 
b and c avoid the 3:1 resonance and be stable at all times cannot be satisfied if Q^ = 21.5 
and k = 0.3. 

However, if the two planets have properties similar to analogous bodies in the Solar 
System, then the system would have been stable, and never in a 3:1 resonance. In Fig. [TO] we 
show the evolution of the 2 planets assuming the inner planet is Neptune-hke {Qb = 10^, Rb = 
3.74i?0), and planet c is a typical terrestrial planet {Qc = 100, k = 0.3, Re = 1.8-R®) (Yoder 



^Publicly available at pittp://janus. astro.umd.edu/HNBody 
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1995). In this case, the 3:1 resonance crossing occurred nearly 10 Gyr ago, and the system 
10 Gyr ago is stable (assuming low amplitude apsidal libration). This scenario invokes the 
simplest assumptions, and is therefore the preferable explanation. We conclude that Gl 581 
c is probably a terrestrial planet with Q ~ 100, but was never inside the 50% cloud cover 
EHZ. 



5. Conclusions 

The detection of a terrestrial planet around a low-mass star is insufficient to determine 
that planet's past and future habitability. The tidal forces between planet and star can 
significantly change the orbits, and hence limit the habitable lifetime. Planets detected in 
the HZ with large eccentricities may be bound for hotter temperatures and ultimately a 
global extinction. On the other hand, planets detected interior to the HZ may have been 
habitable in the past. Gl 581 c was most likely not habitable in the past, but § 4.1 shows 
that if its companion planets were on different orbits, past habitability would have been 
possible. 

For plausible physical and orbital properties of hypothetical terrestrial planets, tides 
may evolve the planets' orbits past the inner edge of the HZ over a timescale comparable to 
the age of the Earth. Therefore in order for planets to be habitable long enough for complex 
life to develop, they must form with eccentricities low enough that tides don't eventually 
make them inhospitable to life. Alternatively planets that form beyond lout niay evolve into 
the habitable zone, but the "flux" of planets out, through the inner edge of the EHZ, is 
larger than that of planets in, through the outer edge, because of the very steep dependence 
of da/dt on a. Therefore, on average and assuming uniform distribution in a, tides tend to 
reduce the total number of habitable planets around M stars in the galaxy. 

Our work has shown that tides may shorten habitable lifetimes of planets that form 
with eccentricities larger than about 0.5. Although the formation mechanism of such planets 
is unknown, the existence of giant planets with such eccentricities suggests even terrestrial- 
sized planets may form with similar values. However, the majority of known exoplanets 
probably formed with eccentricities below 0.5 (see, e.g. Jackson, Greenberg & Barnes 2007), 
thus most habitable lifetimes will not be signiflcantly shortened due to tides. 

This work predicts that the vast majority of terrestrial planets in HZs around low-mass 
stars ( <O.2M0) will be detected on nearly circular orbits because the time to circularize 
the orbit is less than 1 Gyr, c.f. Fig. [31 Should space-based transit missions like COROT 
and Kepler flnd terrestrial planets on non-circular orbits around old (more than a few Gyr), 
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low-mass stars (< 0.35M*) in small orbits (a <0.1 AU), then that system is likely to contain 
additional companions (Mardling & Lin 2004) which pump up eccentricities due to mutual 
interactions. 

The increasing stellar flux on a planet due to orbital decay could impact its atmospheric 
conditions and composition. Previous models of habitable atmospheres have considered 
planets at a single semi-major axis value {i.e. Kasting et al. 1993; Kasting & Catling 2003; 
Segura et a/. 2005; Tinetti, Rashby & Yung 2006; Kiang et a/. 2007). The tidal change in 
the orbit of a habitable planet may affect the relative abundances of biosignatures, such as 
the simultaneous presence of oxygen and methane (Sagan et al. 1993). Future work on the 
evolution of a habitable atmosphere and the identification of biosignatures should consider 
the change in stellar flux due to tidal evolution. 

When a planet is detected such that it was most likely in the HZ in the past, it would 
be interesting to determine if it did, in fact, support life in the past. In the foreseeable 
future, such a determination will have to be made through analysis of disk-averaged spectra 
of planets (Tinetti et al. 2006; Tinetti, Rashby & Yung 2006; Kaltenegger, Traub & Jucks 
2007). Future work on the evolution of habitable atmospheres should explore the possibility 
of detecting signatures of extinct life on planets around M stars. 

Future work should also address the determination of habitable lifetimes in multiple 
planet systems as well as incorporate higher order corrections to the tidal equations. Such 
improvements will require substantial advancement in our understanding of the deforma- 
tions of solid and gaseous bodies due to tides. The Gl 581 systems demonstrates the need 
to develop a theory that describes two tidally damped orbits that also experience signifi- 
cant mutual interactions (note that Mardling [2007] has derived expressions for systems in 
which 1 planet experiences tidal evolution). However, the results of § 4.2 suggest that such 
improvements are unlikely to change our assessment for the Gl 581 system. 

Planets with masses <0.3M,^ are likely too small to support plate tectonics, which 
are thought to be necessary for life (Williams, Kasting & Wade 1997; Raymond, Scalo & 
Meadows 2007). This work suggests some of these smaller planets may be heated internally 
by tides, such as in the Galilean satellites (Peale, Cassen & Reynolds 1979). In some cases, in 
which thl is long enough for life to develop, the tidal working of the body may supply enough 
heat to drive plate tectonics after radiogenic heating has become negligible. Therefore these 
low-mass planets, previously thought to be uninhabitable, may, in fact, be good locations 
for life. Future work should explore this possibility. 

Perhaps the most distressing aspect of this work (from a SETI perspective) is the pre- 
diction that planets can be habitable long enough for complex life to develop, but then that 
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life is extinguished by tides. Yet this work suggests that such a "tidal extinction" may occur 
on some planets around low-mass stars. 

The results presented here may have significant implications in the context of the Gaia 
hypothesis (Lovelock & Margulis 1974), which proposes that biology can change the physical 
conditions on a planet, such as its atmosphere, so as to promote the continuing existence 
of life. In other words, evolution favors organisms that help maintain a habitable environ- 
ment. According to that hypothesis, an inhabited planet may avert a tidal extinction due 
to its organisms developing adaptations that not only change themselves, but also change 
the atmosphere such that the sterilizing effects of increased stellar radiation are mitigated. 
Therefore the Gaia hypothesis predicts that we may observe biosignatures on planets that 
are on orbits interior to the nominal HZ, if they were formerly in the HZ. In that case we 
might need to refine the definition of the HZ into two branches: the "physical" HZ {i.e. 
the definition proposed by Kasting, Whitmire & Reynolds [1993] that applies to planets be- 
fore life gains a foothold), and the "biologically-extended" HZ {i.e. a larger region in which 
evolved life can stave off sterilization). For example, planet Gl 581 c was probably never in 
the physical HZ, but if it once had been, then it might have remained habitable even after 
reaching its current orbit, according to the Gaia hypothesis. The detection of other similar 
new planets may provide observational tests of the Gaia hypothesis. The detection of such a 
planet would be a landmark event, not only in the detection of life beyond the Solar System, 
but because it will reveal fundamental features of life, evolution, and planetary habitability. 

We conclude that tides can move habitable planets into non-habitable orbits and de- 
scribed how this phenomenon is relevant for astrobiological models. Future work into the 
search for habitable worlds, through planet-hunting, atmospheric modeling, geophysical mod- 
eling, etc., should bear in mind the effects of tidal evolution of planets around M stars. 

This work was funded by PG&G grant NNG05GH65G. S.N.R. was supported by an 
appointment to the NASA Postdoctoral Program t the University of Colorado Astrobiology 
Center, administered by Oak Ridge Associated Universities through a contract with NASA. 
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Table 1: Best-Fit Orbital Elements of the Gl 581 Planetary System 
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Figure 1 - Semi-major axes that are habitable as a function of stellar mass and eccentricity 
from Eqs. fl8]- [TTl) . The shading represents different choices of cloud cover for the limits of 
the HZ. Darkest gray assumes no cloud cover, medium gray assumes 50% cloud cover, light 
gray assumes total cloud cover, and white is uninhabitable. 

Figure 2 - Contours of thl of terrestrial planets, that begin at the inner edge of the 0% 
cloud cover EHZ, as a function of stellar mass and initial eccentricity cq. The terrestrial 
planets have a mass of 5 M^ (top left), 1 M^ (top right), 0.5 M® (bottom left), and 0.1 M^ 
(bottom right), and a radius determined by Eq. ([5]). Planets with initial eccentricities below 
0.53, or orbiting stars with a mass greater then 0.35 Mq are habitable for at least 10 Gyr. 
Lower mass planets tend to have longer habitable lifetimes. 

Figure 3 - Evolution of the semi-major axes (solid line) for two hypothetical planets relative 
to EHZ boundaries (shown by the shading, c.f. Fig. [1]). These planets are in orbit about a 
0.15 Mq star. Top: Evolution of a planet with cq = 0.52. This planet's semi-major axis 
decreases until its eccentricity reaches zero after 750 million years, at which point the tidal 
evolution effectively stops. The tidal evolution stopped just before the planet crossed the 
inner edge of the EHZ, and the planet will be habitable indefinitely. Bottom: Evolution of 
a planet with eo = 0.54. Although very similar to the other planet's evolution, this planet's 
larger initial eccentricity results in more torque so that when e reaches zero, the planet is 
just interior to the 50% cloud cover HZ. The habitable lifetime for this planet is 7.5 x 10^ 
years. 

Figure 4 - Top: Orbit-averaged stellar flux received by a tidally evolving planet relative to 
that of the Earth F^ (solid line), see Eq. ([7]), as well as at periastron and apoastron (dotted 
lines). The planet orbits a 0.2 Mq star with initial orbital elements of Oq = 0.0877, and 
Co = 0.84. The tidal lifetime of the planet, 4.5 Gyr in this case, corresponds to the age of the 
Earth, and is shown by the dashed line. Note, however, the extreme difference between flux 
at periastron and apoastron (a factor of 200). Bottom: Ratio of the rotation frequency of the 
tidally evolving planet to its instantaneous mean motion, c.f. Eq. ([3]). Initially the planet 
rotates nearly 10 times faster than it revolves, but after ~ 6 Gyr, the orbit has circularized 
and VL = n. 

Figure 5 - Time since Gl 581 c would have been habitable as a function of its current ec- 
centricity. We set Re = l.SR^, and consider Q values of 100 and 21.5. The Q = 21.5 case 
permits habitability as recently as 3 Gyr ago. 
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Figure 6 - Evolution of Gl 581 c for three different possible radius values. The shading 
represents different definitions of the HZ, c.f. Fig. [1] 

Figure 7 - Evolution of at and e^ for various choices of Qt,, assuming k = 3/2. The evolution 
is negligible unless Qb <10^. 

Figure 8 - Apsidal evolution of the Gl 581 b and c planets, see Table 1. Top: The apsidal 
oscillation is circulation. Bottom: Evolution of Cb (solid line) and Cc (dashed line). The 
orbits are currently near the extrema of their eccentricity values. 

Figure 9 - Evolution of planets b (dashed line), c (solid line) and the instantaneous location 
of b's outer 3:1 mean motion resonance (dotted line) for Qb = 4 x 10^, k = 3/2. In this case 
planet c would have just avoided the 3:1 resonance. The shading corresponds to different 
definitions of the HZ, c.f. Fig. [H 10 Gyr ago the configuration was unstable. 

Figure 10 - Evolution of planets b (short-dashed line), c (solid line) and the instantaneous 
location of b's outer 3:1 mean motion resonance (dotted line) for Qb = 10^, Rb = 3.74i?0, k = 
3/2 and Qc = 100, Re = 1.8i?^, k = 0.3. In our model, this evolution is the most believable, 
and prevents planet c from crossing the 3:1 mean motion resonance and reaching the 50% 
cloud cover habitable zone (medium grey). This evolution permits stable interactions 10 Gyr 
ago. 
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